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Abstract 

There exists a simple rule by which path integrals for the motion of a point particle 
in a flat space can be transformed correctly into those in curved space. This rule 
arose from well-established methods in the theory of plastic deformations, where crys- 
tals with defects are described mathematically by applying nonholonomic coordinate 
transformations to ideal crystals. In the context of time-sliced path integrals, this 
has given rise to a quantum equivalence principle which determines the measure of 
fluctating orbits in spaces with curvature and torsion. The nonholonomic transfor- 
mations are accompanied by a nontrivial Jacobian which in curved spaces produces 
an additional energy proportional to the curvature scalar, thereby canceling an equal 
term found earlier by DeWitt from a naive formulation of Feynman's time-sliced path 
integral in curved space. The importance of this cancelation has been documented 
in various systems (H-atom, particle on the surface of a sphere, spinning top). Here 
we point out its relevance in the process of bosonizing a nonabelian one-dimensional 
quantum field theory, whose fields live in a flat field space. Its bosonized version 
is a quantum-mechanical path integral of a point particle moving in a space with 
constant curvature. The additional term introduced by the Jacobian is crucial for 
the identity between original and bosonized theory. 

A useful bozonization tool is the so-called Hubbard-Stratonovich formula for 
which we find a nonabelian version. 
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I. INTRODUCTION 



In 1956, Bryce DeWitt proposed a path integral formula in curved space using a specific 
generalization of Feynman's time-sliced formula in cartesion coordinates Q. Surprisingly, his 
amplitude turned out to satisfy a Schrodinger equation different from what had previously 
been considered as correct 0: Apart from the Laplace-Beltrami operator for the kinetic 
term, he obtained an extra effective potential proportional to the scalar curvature R. At 
the time of his writing, DeWitt could not think of any argument to outrule the presence of 
such an extra term. 

DeWitt's work has had many successors ||. These employed various time-slicings of the 
action, most popular being postpoint, midpoint, and prepoint prescriptions {§, and added 
to it different correction terms proportional to h 2 to arrive at a Schrodinger equation of their 
personal preference. 

In my opinion, such additional 7z 2 -terms must be rejected since they violate the basic 
principle of Feynman path integrals, according to which a quantum mechanical amplitude 
should be obtainable from a sum over all paths with an amplitude which is the exponential 
gtAu/f^ where A c \ is the purely classical action || along the path. 

The apparent freedom in writing down various path integrals has its counterpart in the 
apparent freedom of setting up a time-evolution operator H from a classical action 

A = J dtLtftf), q^t)) = J dt^g, u (qm u W 

whose Hamiltonian contains products of momenta p^ = dL(g^(£), <^*(i)) / '<9g M and positions 
qt i. 

H(p,q) = ^g^(q)p»P„. (2) 

The metric g^ v {q) describes the geometry of configuration space. If the momenta p^ are 
postulated to satisfy canonical commutation rules with the positions q u , there are many 
different operator orderings corresponding to the same H(p,q). This problem has become 
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known as the operator- ordering problem, and its existence has caused a wide-spread myth 
among theoreticians, that it is basically unsolvable. In fact, many people have expressed 
their belief to the author that different physical systems might have to be quantized with 
different operator orderings. 

Since some years, the author has been fighting this myth. There are many physical 
systems with a Hamiltonian of the form @ for which we know a time-evolution operator 
H whose correctness has never been questioned. The most elementary example is the sym- 
metric spinning top. If the classical Hamiltonian is written as in Eq. (g), with q^ being 
the three Euler angles, and if p^ abd q^ are quantized canonically, there is of course an 
ordering problem. This, however, is due to having chosen the wrong classical variables for 
quantization. Since the system is invariant under rotations but not under translations, only 
the operators associated with the angular momenta Li have good quantum numbers, not 
the generators of translations p^. One must therefore rewrite the classical Hamiltonian (0) 
as [| 

H = (3) 
and impose commutation rules upon the angular momenta Lf. 

[Li,Lj] = ie ijk L k . (4) 

There is no operator-ordering problem in this procedure! 

The same uniqueness holds for any Hamiltonian which is a linear combination of Casimir 
operators and generators of a group of motion in a curved configuration space. These 
observations form the basis of the so-called geometric quantization 0, in which there is no 
source for an extra i?-term. 

Thus we are faced with the problem of finding a construction procedure for path integrals 
in curved spaces which is naturally capable of reproducing these well-established results of 
group quantization. Since path integrals are formulated in phase space in terms of p^ and 
g^-variables which should not be used as a basis for quantization in the operator formulation, 



this seems to be a hard task. Nevertheless, a solution has been found in the form of a simple 
geometric mapping principle. The necessity for finding this solution came from the desire 
to solve the time-sliced path integral of the hydrogen atom, a task which was completed 
in a continumum formulation 17 years ago |J. This solution proceeds by a three-step 
transformation to the path integral of a harmonic oscillator . In the language of ordinary 
quantum mechanics, the three steps proceed as follows: 

First, the Hamiltonian is extended by a dummy forth momentum and written as 

4 V 2 e 2 

Second, a nonholonomic Kustannheimo-Stiefel transformation to coordinates with r = 

2 _ Y^4 /„,u\2 • 



U 



J2fj,=i ( ufl ) is use d to transform H to 



4 u 2 2 

pi — ^ 

This is of the form (|2|) and describes a system in a space with curvature. Only recently it 
was discovered that as a consequence of the nonholonomic nature of the transformation, the 
w^-space carries also torsion |9|-|TT||. 

Classical orbits satisfy energy conservation 

H - E = 0. (7) 

In a third step, the classical equation @ is multiplied by u 2 and becomes 

4 „u2 

V ^- + e 2 - u 2 E = 0. (8) 
^ 8M v ; 

This has a unique operator version describing a harmonic oscillator. 

The intermediate Hamiltonian (|6|) is associated with a unique path integral in a space 
with curvature and torsion, and thus constitutes an important testing ground for any theory. 
If De Witt's construction rules for a path integral in curved space are generalized to such a 
space, one obtains a very complicated Hamilton operator which does not yield the correct 



hydrogen spectrum ]12 



A resolution of this puzzle became possible by the recent discovery of a simple rule 
for correctly transforming Feynman's time-sliced path integral formula from its well-known 
cartesian form corresponding to @ to the spaces with curvature and torsion where the 
dynamics is governed by (||). The rule promises to play a similar fundamental role in 
quantum physics as Einstein's equivalence principle in classical physics, where it fixes the 
form of the equations of motion in curved spaces. It has therefore been named quantum 
equivalence principle (QEP) f|. 

The crucial place where this principle makes a nontrivial statement is in the measure of 
the path integral. The nonholonomic nature of the differential coordinate transformation 
gives rise to an additional term with respect to the naive DeWitt measure, and this cancels 
precisely the bothersome additional term proportional to R in the Schrodinger equation in 
curved space found by DeWitt Q, as well as the many similar additional terms which would 
appear when generalizing De Witt's procedure to spaces with torsion 

It should be mentioned that QEP has drastic consequences even at the classical level if 
the space geometry possesses torsion. As we shall see below, the familiar action principle is 
no longer valid and requires modification: In the presence of torsion, the classical trajectories 
are autoparallels, not geodesies P~5|,|9|]. This surprizing result is most easily illustrated by 



deriving the Euler equations for the motion of a spinning top from an action principle 
formulated within the body-fixed reference frame, where the geometry of the nonholonomic 



coordinates possesses torsion [ |14|j . 



The purpose of this paper is to present another important evidence for the correctness 
of the quantum equivalence principle which arises the context of an exact bosonization of a 
nonabelian fermion model in quantum mechanics. The Hamiltonian of this model is simply 
proportional to square of the total spin of a set of fermions at a point. It is described by a 
set of fermionic harmonic oscillator fields living in a flat field space. When bosonizing this 
model, the fermion fields are replaced by fluctuating angular fields living in a space with 
constant curvature. The associated path integral can be solved 0. The identity between 
initial and bosonized theory give a compelling confirmation for the presence of the nontrivial 



Jacobian generated by the nonholonomic transformation of the path integral measure. 

A similar nonabelian model has, incidentally, been bosonized some 20 years ago ]0| by 
the author in a study of pairing forces in nuclear physics [|TIJ. These forces are described 
by a BCS-like Hamiltonian similar to the one giving rise to the solid-state phenomenon of 
superconductivity. The BCS theory itself was approximately bosonized near the critical 
point almost 40 years ago by Gorkov in his famous derivation of the Ginzburg-Landau 
|17|,[18| . This procedure has been translated into a path integral language almost 20 years 



ago, after developing formalism [|l5j which has since become the prototype for many similar 
enterprises. There exists now a simple theory of collective quantum fields for a wide variety 
of many-body systems, including quarks and gluons . 



The derivation of a Ginzburg-Landau-like theory for superfluid 3 He [f2l3| , |T5|| , and a plas- 



mon description of electron gases [15 were other important applications [^T . 

In superfluid 3 He, the derivation had a novel feature: It was an approximate bosonization 
of a nonabelian system. In order to understand some typical problems arising from the 
nonabelian structure, the author studied in jnj the simple soluble fermion model of nuclear 
pairing forces which he was able to bosonize exactly, ctrrivins* £it el Let grangian of a spinning 
top. However, this bosonization was performed purely formally, without a careful treatment 
of the nonholonomic field transformation whose special properties were unknown at that 
time. The correct result was obtained only by omitting a proper examination of possible 
time slicing corrections. These would have been found to add to the energy an undesirable 
DeWitt type of term proportional to R. 

The recent progress in dealing with nonholonomic field transformations of path integrals 
described in Ref. || enables us to do better. We shall demonstrate that only by performing 
the nonholonomic field transformation according to the new rules provided by the quantum 
equivalence principle does the bosonized theory coincide with the original fermion theory. 

The paper will start in Section II with the bosonization of a rather trivial model, which 
serves to illustrate several essential features of all bosonization procedures. The nonabelian 
model is treated in Section III. 



An important tool for performing abelian bosonizations is the so-called Hubbard- 



Stratonovich transformation formula [22]. Our nonabelian procedure provides us with a 
nonabelian version of this. This formula should be useful for the bosonization of other 
theories, and will be given in Section IV. 

Our results may have consequences for path integral bosonizations of two-dimensional 
nonabelian fermion theories |23|, whose abelian versions were first treated by Coleman, 
Mandelstam, and others |24|. 

Let us first, however, recall the foundations of the quantum equivalence principle. For the 
sake of generality, we shall allow the nonholonomic coordinate transformations to generate 
torsion, just as in the theory of defects, although this general formulation is not required for 
the bosonization to be performed in this paper. 



II. CLASSICAL MOTION OF A MASS POINT IN A SPACE WITH TORSION 

We begin by recalling that Einstein formulated the rules for finding the classical laws of 
motion in a gravitational field on the basis of his famous equivalence principle. He assumed 
the space to be free of torsion since otherwise his geometric priciple was not able to determine 
the classical equations of motion uniquely. Since our nonholonomic mapping principle is not 
beset by this problem, we do not need to rescrict the geometry in this way. The correctness 
of the resulting laws of motion is exemplified by several physical systems with well-known 
experimental properties. Basis for these "experimental verifications" will be the fact that 
classical equations of motion are invariant under nonholonomic coordinate transformations. 
Since it is well known [^,^§] that such transformations introduce curvature and torsion 
into a parameter space, such redescriptions of standard mechanical systems provide us with 
sample systems in general metric-affine spaces. 

To be as specific and as simple as possible, we first formulate the theory for a nonrel- 
ativistic massive point particle in a general metric-affine space. The entire discussion may 
easily be extended to relativistic particles in spacetime. 
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A. Equations of Motion 



Consider the action of the particle along the orbit x(i) in a flat space parametrized with 
rectilinear, Cartesian coordinates: 

A= / dt—{x 1 ) 2 , i = 1,2,3. (9) 
Jt a 2 

It may be transformed to curvilinear coordinates g M , fj, — 1, 2, 3, via some functions 

x* = x\q), (10) 

leading to 

A= / dt—g, u (q)q»q\ (11) 

Jt a ^ 

where 

gM = d,x\q)d v x\q) (12) 

is the induced metric for the curvilinear coordinates. Repeated indices are understood to be 
summed over, as usual. 

The length of the orbit in the flat space is given by 

I = [ tb dty/g fa ,(q)q»q». (13) 

Both the action ( fLl|) and the length fll3|) are invariant under arbitrary reparametrizations of 
space q^ — > q ,fl . 

Einstein's equivalence principle amounts to the postulate that the transformed action 
( p[) describes directly the motion of the particle in the presence of a gravitational field 
caused by other masses. The forces caused by the field are all a result of the geometric 
properties of the metric tensor. 

The equations of motion are obtained by extremizing the action in Eq. ( |TTD with the 
result 



S 



dt(9^q u ) ~ \d^q X q v = + Fj^gV = 0. (14) 

Here 

f \un = 2 (P^9un + d v g\n - d^gxu) (15) 

is the Riemann connection or Christoffel symbol of the first kind. Defining also the Christof- 
fel symbol of the second kind 

VJ* = g^f Xua , (16) 

we can write 

F + f x fq X q" = 0. (17) 

The solutions of these equations are the classical orbits. They coincide with the extrema 
of the length of a curve I in (|i"3|) . Thus, in a curved space, classical orbits are the shortest 
curves, called geodesies. 

The same equations can also be obtained directly by transforming the equation of motion 
from 

a? = (18) 
to curvilinear coordinates g M , which gives 

* = W*' + We^*' = - (19) 

At this place it is useful to employ the so-called basis triads 

dx 1 



*M = 7^7 (20) 



and the reciprocal basis triads 



which satisfy the orthogonality and completeness relations 



- S (21) 



e/eV = 6" v , (22) 
e^e% = 5/'. (23) 

The induced metric can then be written as 

9M = e\{qy u {q)- (24) 

Labeling Cartesian coordinates, upper and lower indices % are the same. The indices //, z/ of 
the curvilinear coordinates, on the other hand, can be lowered only by contraction with the 
metric g^ v or raised with the inverse metric g^ u = (g^u) -1 . Using the basis triads, Eq. Q19[) 
can be rewritten as 

j f (e\qn = e\r + e\,„r? = 0, (25) 

or as 

<f + efe\, x q K q x = 0. (26) 

The subscript A separated by a comma denotes the partial derivative d\ = d/dq x , i.e., 
f x = d\f. The quantity in front of q K q x is called the affine connection: 

Tx^ = efe\,\- (27) 
Due to (|2"2l), it can also be written as 

r A / = -e\e t \x- (28) 
Thus we arrive at the transformed flat-space equation of motion 

<f + T KX »q K q x = 0. (29) 

The solutions of this equation are called the straightest lines or autoparallels . 

If the coordinate transformation functions x l (q) are smooth and single- valued, they are 
integrable, i.e., their derivatives commute as required by Schwarz's integrability condition 
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(d x d K -d K d x )x i (q) = o. (30) 

Then the triads satisfy the identity 

<a = 4 >K , (31) 

implying that the connection is symmetric in the lower indices. In this case it coincides 
with the Riemann connection, the Christoffel symbol F^. This follows immediately after 
inserting g^q) = e\(q)e\(q) into (^) and working out all derivatives using ([!!]). Thus, 
for a space with curvilinear coordinates q^ which can be reached by an integrable coordinate 
transformation from a flat space, the autoparallels coincide with the geodesies. 

B. Nonholonomic Mapping to Spaces with Torsion 

It is possible to map the x-space locally into a g-space via an infinitesimal transformation 

dx* = e\(q)dqv, (32) 

with coefficient functions e\(q) which are not integrable in the sense of Eq. (|30"D, i.e., 

d^^q) - cUV(<?) ± 0. (33) 

Such a mapping will be called nonholonomic. It does not lead to a single- valued function 
x l (q). Nevertheless, we shall write (p3|) in analogy to (p0|) as 

(d x d K -d K d x )x\q)^0, (34) 

since this equation involves only the differential dx l . Our departure from mathematical 
conventions will not cause any problems. 

From Eq. (^) we see that the image space of a nonholonomic mapping carries torsion. 
The connection T Xk ^ = e^e 1 K>x has a nonzero antisymmetric part, called the torsion tensor-^ 



Our notation for the geometric quantities in spaces with curvature and torsion is the same as in 
J. A. Schouten, Ricci Calculus, Springer, Berlin, 1954. 
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s A / = ^(r A /-r K /). (35) 

In contrast to r A(t M , the antisymmetric part S Xk ^ is a proper tensor under general coordinate 
transformations. The contracted tensor 

S, = S, x x (36) 

transforms like a vector, whereas the contracted connection r„ = Y ^ does not. Even 



^ ic m o.^ o foncnr Tiro cVioll ^Voolir lrvtiroT* QTirl Tdica 



though Y^ V X is not a tensor, we shall freely lower and raise its indices using contractions 



with the metric or the inverse metric, respectively: T^ u x = g^ K T KU x , T fl uX = g UK T^ K x , 
T^x = 9\k^hv K ■ The same thing will be done with f^j/. 

In the presence of torsion, the connection is no longer equal to the Christoffel symbol. 
In fact, by rewriting T^ u \ = e^d^v trivially as 



r^A = ^ { e iA<9 M eV + d^e iX e l u + e iiX d v e\ + d v e iiX e\ - e ifl d x e l u - <9 A e^eV} 
+ \ { [ea^eV - eu^e'J - [e iM ^e l A - e^<9 A eV] + [e iu d x e\ - e it ,<9 M e l A ] } 



2 

2 ^ 

and using e\(q)e l u (q) = g^io), we find the decomposition 



V = V + K, v \ (37) 
where the combination of torsion tensors 

K^ uX = — + S^i, (38) 
is called the contortion tensor. It is antisymmetric in the last two indices so that 

In the presence of torsion, the shortest and straightest lines are no longer equal. Since 
the two types of lines play geometrically an equally favored role, the question arises as to 
which of them describes the correct classical particle orbits. The answer will be given at the 
end of this section. 
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The main effect of matter in Einstein's theory of gravitation manifests itself in the viola- 
tion of the integrability condition for the derivative of the coordinate transformation x l (q), 
namely, 

(d^d l/ -d v d fl )d x x i (q)^0. (40) 

A transformation for which x l (q) itself is integrable, while the first derivatives d ll x l {q) = e l ^ 
are not, carries a flat-space region into a purely curved one. The quantity which records the 
nonintegrability is the Cartan curvature tensor 

R,uX K = e^id.d, - cUge\. (41) 

Working out the derivatives using (^7]) we see that R^ uX K can be written as a covariant curl 
of the connection, 

R^v\ K = dfjT u x K — d v v^ — [r M , t u ] x k . (42) 

In the last term we have used a matrix notation for the connection. The tensor components 
T ^ are viewed as matrix elements (r ^) X K , so that we can use the matrix commutator 

[r M , Tv}^ = (r^r^ — r l/ r A1 )A fv = T^ x a v u ^ — y^y^. (43) 

Einstein's original theory of gravity assumes the absence of torsion. The space proper- 
ties are completely specified by the Riemann curvature tensor formed from the Riemann 
connection (the Christoffel symbol) 

R^x = d,T^ - djY^ - [f „, f „] A ". (44) 

The relation between the two curvature tensors is 

Rfiv\ K = R^x + DpK vX « - D V K^ - [Kp, K V \ X K . (45) 

In the last term, the K^'s are viewed as matrices (K^) X K . The symbols denote the 
covariant derivatives formed with the Christoffel symbol. Covariant derivatives act like 
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ordinary derivatives if they are applied to a scalar field. When applied to a vector field, they 
act as follows: 

L>X = + T^v\ (46) 

The effect upon a tensor field is the generalization of this; every index receives a correspond- 
ing additive T contribution. 

In the presence of torsion, there exists another covariant derivative formed with the affinc 
connection V ^ rather than the Christoffel symbol which acts upon a vector field as 

DfjV v = d^v v - T^vx, 

D^ = d^ + T,x u v\ (47) 

This will be of use later. 

From either of the two curvature tensors, R iiV x K and R^ v x K -> one can f° rm the once- 
contracted tensors of rank 2, the Ricci tensor 

Ru\ = R^ux^ 1 , (48) 

and the curvature scalar 

R = 9 uX RuX- (49) 

The celebrated Einstein equation for the gravitational field postulates that the tensor 

Gp,, = R^u — -QfiuR, (50) 

the so-called Einstein tensor, is proportional to the symmetric energy-momentum tensor of 
all matter fields. This postulate was made only for spaces with no torsion, in which case 
R^y = R^ u and R^, are both symmetric. As mentioned before, it is not yet clear 
how Einstein's field equations should be generalized in the presence of torsion since the 
experimental consequences are as yet too small to be observed. In this paper, we are not 
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concerned with the generation of curvature and torsion but only with their consequences 
upon the motion of point particles. 

Two nonholonomic sample mappings producing curvature and torsion are shown in Fig. 
1. They are used in the theory of defects to produce a crystal with a single dislocation or 
disclination, respectively. Readers not familiar with this subject are advised to consult the 
Refs. |25|,|26| and the previous literature on this subject quoted therein. 



Consider first the upper example in which a dislocation is generated, characterized by a 
missing or additional layer of atoms (see Fig. 10.1). In two dimensions, it may be described 
differentially by the transformation 

{dq 1 for i — 1, 

(51) 
dq 2 + ed^ l (f)(q)dq' J - for i = 2, 

with the multi-valued function 

(f)(q) = arctan(g 2 /? 1 )- (52) 
The triads reduce to dyads, with the components 

e\ = 5\ + ed^{q) , (53) 
and the torsion tensor has the components 

e\S^ v x = 0, e\S, p x = e -(d,d v - d v d^. (54) 

If we differentiate ( |5*2"D formally, we find (d ll d u — d u d^)<j) = 0. This, however, is incorrect at 



the origin. Using Stokes' theorem we see that 

/ d 2 q(d 1 d 2 - d 2 d 1 )^> = j> dq^d^ = I d<\> = 2tt (55) 

for any closed circuit around the origin, implying that there is a 5-function singularity at 
the origin with 
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e\S 12 X = ^S {2) (q). (56) 

By a linear superposition of such mappings we can generate an arbitrary torsion in the q- 
space. The mapping introduces no curvature. When encircling a dislocation along a closed 
path C, its counter image C in the ideal crystal does not form a closed path. The closure 
failure is called the Burgers vector 

V = ^dx i = jd<fe\. (57) 

It specifies the direction and thickness of the layer of additional atoms. With the help of 
Stokes' theorem, it is seen to measure the torsion contained in any surface S spanned by C: 

V = I d 2 s^d^\ = / d 2 s^S^ u \ (58) 

J S J s 

where d 2 s fiu = —d 2 s UfJ- is the projection of an oriented infinitesimal area element onto the 
plane \w. The above example has the Burgers vector 

6*=(0 |C ). (59) 

A corresponding closure failure appears when mapping a closed contour C in the ideal 
crystal into a crystal containing a dislocation. This defines a Burgers vector: 

V i =<f d<f=<f dzV- (60) 
By Stokes' theorem, this becomes a surface integral 



= -^ S Ve 3 V, (61) 



the last step following from 

The second example is the nonholonomic mapping in the lower part of Fig. 1 generating 
a disclination which corresponds to an entire section of angle a missing in an ideal atomic 
array. For an infinitesimal angel a, this may be described, in two dimensions, by the 
differential mapping 
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x' = 5\{q^ + ne^q^{q)l (62) 

with the multi- valued function (|52]). The symbol e^ v denotes the antisymmetric Levi-Civita 
tensor. The transformed metric 

2Vt 

9»u = V - —e^ x e\q x q K - (63) 
<fq a 

is single-valued and has commuting derivatives. The torsion tensor vanishes since (8182 — 
d2di)x 1,2 is proportional to g 2 ' 1 ^ 2 -'^) = 0. The local rotation field uj(q) = \{d\x 2 — cfax 1 ), on 
the other hand, is equal to the multi- valued function — 00(g), thus having the noncommuting 
derivatives: 

{d 1 d 2 - d 2 d 1 )u;(q) = -2nQ5^(q). (64) 

To lowest order in Q, this determines the curvature tensor, which in two dimensions posses 
only one independent component, for instance Ruu- Using the fact that g^ v has commuting 
derivatives, Ryin can be written as 

#1212 = (<9i<9 2 - d 2 d 1 )u(q). (65) 



C. New Equivalence Principle 

In classical mechanics, many dynamical problems are solved with the help of nonholo- 
nomic transformations. Equations of motion are differential equations which remain valid if 
transformed differentially to new coordinates, even if the transformation is not integrable in 
the Schwarz sense. Thus we postulate that the correct equation of motion of point particles 
in a space with curvature and torsion are the images of the equation of motion in a flat 
space. The equations (|29|) for the autoparallels yield therefore the correct trajectories of 
spinless point particles in a space with curvature and torsion. 

This postulate is based on our knowledge of the motion of many physical systems. Im- 
portant examples are the Coulomb system 0, and the spinning top described with nonholo- 
nomic coordinates within the body-fixed reference system WM. Thus the postulate has a 



good chance of being true, and will henceforth be referred to as a new equivalence principle. 
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D. Classical Action Principle for Spaces 
with Curvature and Torsion 



Before setting up a path integral for the time evolution amplitude we must find an action 
principle for the classical motion of a spinless point particle in a space with curvature and 
torsion, i.e., the movement along autoparallel trajectories. This is a nontrivial task since 
autoparallels must emerge as the extremals of an action ([IT]) involving only the metric tensor 
g^y. The action is independent of the torsion and carries only information on the Riemann 
part of the space geometry. Torsion can therefore enter the equations of motion only via 
some novel feature of the variation procedure. Since we know how to perform variations of an 
action in the euclidean x-space, we deduce the correct procedure in the general metric-affine 
space by transferring the variations Sx l (t) under the nonholonomic mapping 

i? = eS(q)x l (66) 

into the g M -space. Their images are quite different from ordinary variations as illustrated in 
Fig. |X|(a). The variations of the Cartesian coordinates 5x l (t) are done at fixed end points of 
the paths. Thus they form closed paths in the x-space. Their images, however, lie in a space 
with defects and thus possess a closure failure indicating the amount of torsion introduced 
by the mapping. This property will be emphasized by writing the images dq^(t) and calling 
them nonholonomic variations. 

Let us calculate them explicitly. The paths in the two spaces are related by the integral 
equation 

<f W = q^ta) + f dt'ef(q(t'))x\f). (67) 

Jt a 

For two neighboring paths in x-space differing from each other by a variation 5x l (t), Eq. fl5"?D 
determines the nonholonomic variation Sq^lt): 

dq»(t) = t d£d[ei"{q{l!))x\l!)\. (68) 

Jta 

A comparison with ( |5"6j ) shows that the variations <5g M and the time derivative of q^ are 
independent of each other 
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just as for ordinary variations 8x l . 

Let us introduce an auxiliary holonomic variations in g-space: 

5q» = ef{q)8x\ 

In contrast to Sq^(t), these vanish at the endpoints, 

5q(t a ) = 5q(t b ) = 0, 

i.e., they form closed paths with the unvaried orbits. 
Using (|70|) we derive from fl68D the relation 

j t W(t) = de^{q{t))x\t) + enq{t))dx\t) 

= WW))**® +ef{q{t))j t [e\{tW{t)]. 



After inserting 



deS(q) = -IV*<fgV, M = T x ^q x e\, 



this becomes 



^M(t) = -r Ai /cfgV + T^q x 5q v + Uq». 
at at 

It is useful to introduce the difference between the nonholonomic variation Sq^ 
auxiliary holonomic variation Sq^: 

W = Sq^ 1 - 5qV. 

Then we can rewrite ([73]) as a first-order differential equation for Sb* 1 : 

^-W = -T x ^Sb x q" + 2S x ^q X 8q v . 
at 

After introducing the matrices 
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(77) 



and 



W v (t)=2S x ^(q{t))q\t), 



(78) 



equation (|76|) can be written as a vector differential equation: 

d 



This is solved by 



dt 



db= -Gdb+E(t) 5q v (t). 



(79) 



with the matrix 



db(t) = / dt'U(t,t') E(t') 6q(t'), 

Jta 



U(t,t') =Texp [- f dt"G(t 
Jt> 



(80) 



(81) 



In the absence of torsion, E(t) vanishes identically and Sb(t) = 0, and the variations dq^(t) 
coincide with the holonomic 5q^(t) [see Fig. 0(b)]. In a space with torsion, the variations 
dq^(t) and 5q^{t) are different from each other [see Fig. 0(c)]. 

Under an arbitrary nonholonomic variation dq^{t) = 5q^ + <Jfe M , the action changes by 



SA = M ^ dt (g^fdqv + ^d,g x Jq»q x q K 



(82) 



After a partial integration of the <5g-term we use (|7T|), (|69"D, and the identity d il g vX = T^ x + 
r^Ai/; which follows directly form the definitions g^ u = e^eV und = e i A 9 A1 eV, and 



obtain 



dA = M dt 

>t a 



r £" 9flv (q v + f XK u q X q K ) W + (g^j^ + T^ x J¥ l q x q K 



(83) 



To derive the equation of motion we first vary the action in a space without torsion. 
Then db^(t) = 0, and we obtain 



SA = 5A= -M / dtg^if + Y Xk u q q KS )q u . 



(84) 
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Thus, the action principle SA = produces the equation for the geodesies (|T7D, which are 
the correct particle trajectories in the absence of torsion. 

In the presence of torsion, Sb^ ^ 0, and the equation of motion receives a contribution 
from the second parentheses in (^). After inserting Q7E|), the nonlocal terms proportional 
to Sb^ cancel and the total nonholonomic variation of the action becomes 

dA = -M f" dbg. 

Jta 

= -M f b dt 9fU , (f + r XK »q x q K ) (85) 

The second line follows from the first after using the identity f \ K U = T{\ K } U + 2S V '{a k }- The 
curly brackets indicate the symmetrization of the enclosed indices. Setting dA = gives the 
autoparallels (53) as the equations of motions, which is what we wanted to show. 

III. ALTERNATIVE FORMULATION OF ACTION PRINCIPLE 

WITH TORSION 

The above variational treatment of the action is still somewhat complicated and calls for 
an simpler procedure which we are now going to present. f\ 

Let us vary the paths q^{t) in the usual holonomic way, i.e., with fixed endpoints, and 
consider the associated variations 5x' L = e l ^{q)5q^ of the cartesian coordinates. Taking their 
time derivative d t = d/dt we find 

dtdx' = e\(q)d t 5q X + d ll e\(q)q»5q\ (86) 

On the other hand, we may write the relation (^) in the form d t x l = e l fJj (q)d t q fJ- and vary 
this to yield 

5d t x l = e\(q)5q x + d,e\(q)q x 5q». (87) 
Using now the fact that time derivatives 5 and variations dt commute for cartesian paths, 



2 See H. Kleinert und A. Pelster, FU-Berlin preprint, May 1996. 



<f + (f a/ + 2S\ K ) q X q K 5q 
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5d t x i -dtSx 1 = 0, (88) 
we deduce from (|8q) and (87) that this is no longer true in the presence of torsion, where 



6d t q x -d t 6q x = 2S^(q)r6q». (89) 

In other words, the variations of the velocities q^{t) no longer coincide with the time deriva- 
tives of the variations of q^(t). 

This failure to to commute is responsible for shifting the trajectory from geodesies to 
autoparallels. Indeed, let us vary an action 

t b 

A = J dtL(q x (t),q x (t)) (90) 



by 5q x (t) and impose ([39"1) , we find 

s H d iw^ + w^ +2S -w (91) 

After a partial integration of the second term using the vanishing 5q x (t) at the endpoints, 
we obtain the Euler-Lagrange equation 

9L d9 \ = 2S x ^ (92) 



dq x dtdq x dq v 

This differs from the standard Euler-Lagrange equation by an additional contribution due 
to the torsion tensor. For the action (|TTD we thus obtain the equation of motion 

M 



9 XK (d,g VK - \ 8 K g, v ) + 2S\ u ]q»cf = 0, (93) 



which is once more Eq. (29) for autoparallels. 



IV. PATH INTEGRAL IN SPACES WITH CURVATURE AND TORSION 

We now turn to the quantum mechanics of a point particle in a general metric-affine 
space. Proceeding in analogy with the earlier treatment in spherical coordinates, we first 
consider the path integral in a flat space with Cartesian coordinates 
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fxtlx't") 



-D 



N 

n 



N+l 



n ^o(Ax n ), 



71=1 



'27rieh/M n=i 

where ^Q(Ax n ) is an abbreviation for the short-time amplitude 



(94) 



i^(Ax n ) = (x n | exp ( ) |x n _!) 



1 



-D 



exp 



i M (Ax r 



2 



(95) 



'2meh/M 

with Ax„ = x n — x n _i, x = xjv+i, x' = Xq. A possible external potential has been omitted 
since this would contribute in an additive way, uninfluenced by the space geometry. 

Our basic postulate is that the path integral in a general metric-affine space should be 
obtained by an appropriate nonholonomic transformation of the amplitude fl94"| ) to a space 
with curvature and torsion. 



A. Nonholonomic Transformation of the Action 

The short-time action contains the square distance (Ax n ) 2 which we have to transform 
to g-space. For an infinitesimal coordinate difference Ax n « dx n , the square distance is 
obviously given by (dx.) 2 = g^dq^dq" . For a finite Ax n , however, it is well known that 
we must expand (Ax n ) 2 up to the fourth order in Ag n M = q n ^ — g n _i M to find all terms 
contributing to the relevant order e. 

It is important to realize that with the mapping from dx 1 to dq^ not being holonomic, 
the finite quantity Ag M is not uniquely determined by Ax 1 . A unique relation can only be 
obtained by integrating the functional relation ( |B7| ) along a specific path. The preferred 
path is the classical orbit, i.e., the autoparallel in the g-space. It is characterized by being 
the image of a straight line in the x-space. There x l (t) =const and the orbit has the linear 
time dependence 

Ax\t) = i i (t )At, (96) 

where the time to can lie anywhere on the t-axis. Let us choose for i the final time in each 
interval (t n ,t n -i). At that time, x l n = x l {t n ) is related to q£ = g M (t n ) by 
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4 = ^Mtt- 



(97) 



It is easy to express qf^ in terms of Aq% = q£ — q^-i along the classical orbit. First we expand 
g M (t n _i) into a Taylor series around t n . Dropping the time arguments, for brevity, we have 



e 2 e 3 a 
Aq = q X - q' X = eq X - -q X + -q + ... , 



(98) 



where e = t n — t n -\ and q x , q x , . . . are the time derivatives at the final time t n . An expansion 
of this type is referred to as a postpoint expansion. Due to the arbitrariness of the choice 
of the time to in Eq. (0), the expansion can be performed around any other point just as 
well, such as and t n = (t n + t n _i)/2, giving rise to the so-called prepoint or midpoint 
expansions of Aq. 

Now, the term q x in ( |98|) is given by the equation of motion (p9|) for the autoparallel 



(99) 



A further time derivative determines 



if = -(d CT rv A - 2r M /r {(TT} A )^^g CT 



(100) 



Inserting these expressions into (|98| ) and inverting the expansion, we obtain q x at the final 
time t n expanded in powers of Aq. Using (|96| ) and (|9~?D we arrive at the mapping of the 
finite coordinate differences: 



Ax 1 = e\q x At 



(101) 



Aq x -^T, u x Aq»Aq» + j ] {d a T, v x + T^T {aT} X )Aq»Aq»Aq° + . . . 



where e l \ and are evaluated at the postpoint. Inserting this into the short-time ampli- 



tude (|95|) , we obtain 



fq(Ax) = (x| exp ( --eH ) |x - Ax) : 



with the short-time postpoint action 



1 



-D 



exp 



2irieh/M 



-Al{q,q-Aq) 



(102) 
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M ■ o M 
A%(q,q- Ag) = ^(Ax 1 ) 2 = e-g^? 

= j^Ag" Atf - r^ x Aq»Aq»Aq x (103) 
+ \^9„r{d K Yxu T + r Xu s T {K5} T ) + ^T^] Ag M Ag I/ Ag A Ag K + ... }. 

Separating the affine connection into Christoffel symbol and torsion, this can also be written 

as 

Mr 



AUQ, q~Aq) = — {^A^Ag" - T^Aq^Aq 



+ 



Aq' 1 Aq u Aq x Aq K 



1 



(104) 



Note that the right-hand side contains only quantities intrinsic to the g-space. For the 
systems treated there (which all lived in a euclidean space parametrized with curvilinear 
coordinates), the present intrinsic result reduces to the previous one. 

At this point we observe that the final short-time action ( |103| ) could also have been 
introduced without any reference to the flat reference coordinates x 1 . Indeed, the same 
action is obtained by evaluating the continuous action (|TTj) for the small time interval At = e 
along the classical orbit between the points g n _i and q n . Due to the equations of motion 



(f29|), the Lagrangian 



%!) = Yff^('))W(*) 



(105) 



is independent of time (this is true for autoparallels as well as geodesies). The short-time 
action 

M ft 



A^q,q') 



dtg,Mt)mt)<i u (t) 



2 Jt-e 

can therefore be written in either of the three forms 

A e = —eg^{q)q"q V = ye^OzW = —e 9liV {q)<f<f, 



(106) 



(107) 



where g M , g' M , g M are the coordinates at the final time t n , the initial time £ n _i, and the average 



time (t n +t„_i)/2, respectively. The first expression obviously coincides with (|107|) . The 
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others can be used as a starting point for deriving equivalent prepoint or midpoint actions. 
The prepoint action arises from the postpoint one At, by exchanging Aq by — Aq and the 
postpoint coefficients by the prepoint ones. The midpoint action has the most simple-looking 
appearance: 

A e (q+^,q~^)= (108) 



M 
Ye 



g^^^Q U + ^r(dxT^ T + T^T {x5} T )Aq^Aq u Aq x Aq K + 



where the affine connection can be evaluated at any point in the interval {t n -i,t n ). The 
precise position is irrelevant to the amplitude producing only changes beyond the relevant 
order epsilon. 

We have found the postpoint action most useful since it gives ready access to the time 
evolution of amplitudes, as will be seen below. The prepoint action is completely equivalent 
to it and useful if one wants to describe the time evolution backwards. Some authors favor the 
midpoint action because of its symmetry and intimate relation to an ordering prescription in 
operator quantum mechanics which was advocated by H. Weyl. This prescription is, however, 
only of historic interest since it does not lead to the correct physics. In the following, the 
action A e without subscript will always denote the preferred postpoint expression (|103 ): 



A e = A%{q,q- Aq). (109) 



B. The Measure of Path Integration 

We now turn to the integration measure in the Cartesian path integral (fM|) 

1 

, =D IT ^ Xn - 

^2meh/M n=i 

This has to be transformed to the general metric-affine space. We imagine evaluating the 
path integral starting out from the latest time and performing successively the integrations 
over xn,xn-i, • • • , i-e., in each short-time amplitude we integrate over the earlier position 
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coordinate, the prepoint coordinate. For the purpose of this discussion, we relabel the 
product n^Li dPx l n by Y^,=i dx % n _ x , so that the integration in each time slice with 
n = N + 1, N, . . . runs over dx l n _ l . 

In a flat space parametrized with curvilinear coordinates, the transformation of the 
integrals over d D x l n _ 1 into those over d D q^ f _ 1 is obvious: 



N+l r N+l . 

II / dD <-i = II / d°qU det ej,^ 

n=2 J n=2 ^ J 



The determinant of e l ^ is the square root of the determinant of the metric g^: 



(110) 



det(eV) = Jdet g^(q) = J g{q) 



111) 



and the measure may be rewritten as 

N+l N+l r r -I 

II / d°<-i = II / dD( ln-i yfgijb^) ■ (112) 

n=2 J n=2 LJ J 

This expression is not directly applicable. When trying to do the d D (^^-integrations suc- 
cessively, starting from the final integration over dq 1 ^, the integration variable q n -i appears 
for each n in the argument of det e^(g n _i) or g^^n-i)- To make this g n _i-dependence 
explicit, we expand in the measure ( |110| ) e^(g n _i) = e\(q n — Aq n ) around the postpoint q n 
into powers of Aq n . This gives 

dx* = e;(g - Aq)dq» = e^dq» - e\^Aq v + \e\ )VX dq» Aq" 'Aq x + . . . , (113) 

omitting, as before, the subscripts of q n and Aq n . Thus the Jacobian of the coordinate 
transformation from dx l to dq 11 is 



Jo = det( e g det [5% - efe\, v Aq v + -e.V^A^ Ag A 
giving the relation between the infinitesimal integration volumes d D x % and d D q fl : 

N+l N+l , s 

II d D x\ l _ 1 = 1] { d D q^ l J 0n \. 

71=2 71=2 

The well-known expansion formula 



(114) 



;ii5) 
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det(l + B) = exptrlog(l + B) = exptr(5 - B 2 /2 + B 3 /3 



(116) 



allows us now to rewrite Jn as 



J = det(e* K )exp ( -A\ j . 



(117) 



with the determinant det(e^) = y g(q) evaluated at the postpoint. This equation defines an 
effective action associated with the Jacobian, for which we obtain the expansion 



-A% = -e.V^Ag" + - 



Aq u Aq x + . . . . 



(118) 



To express this in terms of the affine connection, we use (p7|) and derive the relations 



-e- \ = ~e CT e* e- e j = -T a I\ 

^c-i!/,^ ft, A 4 3 U K >X 4 iiv i 1 After 

T^Gi/iC f,A/t 7j9[iT [9k ifii 6 v,a) G G v,X^ ,k\ 



A;/ Tkct )• 



(119) 



(120) 



With these, the Jacobian action becomes 



h Aj ° 



-T^AqV + -d^^Aq'Aq* + ... . 



;i2i) 



The same result would, of course, be obtained by writing the Jacobian in accordance with 
©as 



Jo = yg(q- &q), 



(122) 



which leads to the alternative formula for the Jacobian action 



ex P ( r-4!/o 



g(q - Aq) 



g{q) 



(123) 



An expansion in powers of Aq gives 

1 



1 



9 m v^)^^+- • 



(124) 
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Using the formula 



1 „ _ 1 



d^V9 = ^ <T<W = V, (125) 



this becomes 

exp ( ) = i - r M /A^ + ±(^f , A A +f M /r vX x )£*f&<f + ... 



(126) 



so that 



= -VV + \duF uX X Aq^Aq v + . . . . (127) 
In a space without torsion where T x u = T^ u x , the Jacobian actions ( 121|) and ( |127| ) are 



trivially equal to each other. But the equality holds also in the presence of torsion. Indeed, 



when inserting the decomposition (|3T|), T^ u x = T x + K tlv x : into ( |121| ), the contortion tensor 
drops out since it is antisymmetric in the last two indices and these are contracted in both 
expressions. 

In terms of A € Jon , we can rewrite the transformed measure (|110|) in the more useful form 



N+1 . N+1 



II / d ° x n-i = II ( / dD( in-i det [e;(g n )] exp (U% n )\. (128) 

n=2 J n=2 U " VAi / J 

In a flat space parametrized in terms of curvilinear coordinates, the two sides of (|1 10|) 
and (|1 28|) are related by an ordinary coordinate transformation, and the right-hand side 



gives the correct measure for a time-sliced path integral. In a general metric-affine space, 
however, this is no longer true. Since the mapping dx l — > dq 11 is nonholonomic, there are 
in principle infinitely many ways of transforming the path integral measure from Cartesian 
coordinates to a noneuclidean space. Among these, there exists a preferred mapping which 
leads to the correct quantum-mechanical amplitude in all known physical systems. It is this 
mapping which led to the correct solution of the path integral of the hydrogen atom || . 
The clue for finding the correct mapping is offered by an unesthetic feature of Eq. ( |113| ) : 



The expansion contains both differentials dq^ and differences Aq^. This is somehow incon- 
sistent. When time-slicing the path integral, the differentials dq^ in the action are increased 
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to finite differences Ag M . Consequently, the differentials in the measure should also become 
differences. A relation such as ( |113| ) containing simultaneously differences and differentials 
should not occur. 

It is easy to achieve this goal by changing the starting point of the nonholonomic mapping 
and rewriting the initial flat space path integral 



as 



(xtlxV) 



1 



N 

2nieh/M n=i 



d D Ax r 



N+l 



n ^(Ax n ). 



(129) 



n=l 



Note that since Q n are Cartesian coordinates, the measures of integration in the time-sliced 
expressions (0) and Q129Q are certainly identical: 

N „ N+l 



IJ / d D x n = IJ / d D Ax r , 

n=l n=2 



(130) 



Their images under a nonholonomic mapping, however, are different so that the initial form 
of the time-sliced path integral is a matter of choice. The initial form ( |129| ) has the obvious 
advantage that the integration variables are precisely the quantities Ax l n which occur in the 
short-time amplitude K^{Ax n ). 

Under a nonholonomic transformation, the right-hand side of Eq. ( |130|) leads to the 
integral measure in a general metric-affine space 

N+l . N + l r 



n / d D Ax n - n / d ° A( in Jn 

n=2 J n=2 <- J 



(131) 



with the Jacobian following from ( |101| ) (omitting n 
d(Ax) 

1 



J= 



d(Aq) 

--det(e l K ) det 



(132) 



5/-r w A Ag^ + -(a (T V + r {At /r {T | (T}} A )Ag-Ag CT + . 



In a space with curvature and torsion, the measure on the right-hand side of ( |131|) replaces 



the flat-space measure on the right-hand side of (|112 ). The curly double brackets around 
the indices u, k, a, fi indicate a symmetrization in r and a followed by a symmetrization in 
/i, u, and a. With the help of formula (|116|) we now calculate the Jacobian action 
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+- 



1 



(133) 



{uk} l{a M } 



A<f Ag^ + . . . 



This expression differs from the earlier Jacobian action ( |121| ) by the symmetrization symbols. 
Dropping them, the two expressions coincide. This is allowed if g M are curvilinear coordinates 
in a flat space. Since then the transformation functions x l (q) and their first derivatives 
d fl x 1 (q) are integrable and possess commuting derivatives, the two Jacobian actions (|121[) 
and (|133|) are identical. 

There is a further good reason for choosing ( |130| ) as a starting point for the nonholonomic 
transformation of the measure. According to Huygens' principle of wave optics, each point 
of a wave front is a center of a new spherical wave propagating from that point. Therefore, 
in a time-sliced path integral, the differences Ax l n play a more fundamental role than the 
coordinates themselves. Intimately related to this is the observation that in the canonical 
form, a short-time piece of the action reads 

dp. 



2nh 



exp 



i_ , _ , _ ip 2 n 

h Pn{X n X n - X ) 2Mh t 



(134) 



Each momentum is associated with a coordinate difference Ax n = x n — x n -\. Thus, we 
should expect the spatial integrations conjugate to p n to run over the coordinate differences 
Ax n = x n — x n _i rather than the coordinates x n themselves, which makes the important 
difference in the subsequent nonholonomic coordinate transformation. 

We are thus led to postulate the following time-sliced path integral in g-space: 

1 



(g|exp 



-j-(t-t')H 



W) 



N+l 

i, n 

2-nihe/M n=i 



d D Aq„ 



9(<ln) 



-D 



2irieh/M 



x exp 



11 n=l 



(135) 



where the integrals over Aq n may be performed successively from n = N down to n — 1. 

Let us emphasize that this expression has not been derived from the flat space path 
integral. It is the result of a specific new quantum equivalence principle which rules how a 
flat space path integral behaves under nonholonomic coordinate transformations. 
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It is useful to reexpress our result in a different form which clarifies best the relation 
with the naively expected measure of path integration ( |112j ), the product of integrals 

N „ N 



II Jd D x n =U Jd D q n J^h) 

n=l n=l 



(136) 



The measure in (|135|) can be expressed in terms of (|136|) as 



N+l 

n 

n=2 



d D Aq n Jg(q n 



N f 

n d D qn ^)e- iA y n 

n=l LJ 



The corresponding expression for the entire time-sliced path integral (|135|) in the metric- 
affine space reads 



(g|exp 



-^(t-t')H \q') 



-D 



N 

n 



d D q n 



x exp 



2mhe/M *=i 

N+l 



-D 



2-Kihe/M 



n=l 



(137) 



where AAj is the difference between the correct and the wrong Jacobian actions in Eqs. ( |121| ) 
and P3 ): 



AA e j = A e j - A%. 



(138) 



In the absence of torsion where T 



{/lis} 



r A 



this simplifies to 



(139) 



where R^ u is the Ricci tensor associated with the Riemann curvature tensor, i.e., the con- 



traction (48) of the Riemann curvature tensor associated with the Christoffel symbol V 



Being quadratic in Aq, the effect of the additional action can easily be evaluated per- 
turbatively using the methods explained in Chapter 8, according to which Aq^Aq u may be 
replaced by its lowest order expectation 

(A<f A<f) = ie%g^{q)/M. 
Then AAj yields the additional effective potential 
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cff 



6M 



R, 



(140) 



where R is the Riemann curvature scalar. By including this potential in the action, the path 
integral in a curved space can be written down in the naive form ( |136| ) as follows: 



(g|exp 



1 



N 

n 

2mhe/M n=i 

N+l 



-D 



d D q n 



^/2meh/M D 



x exp 



rE(^ £ + ^ eff ) 



71=1 



;i4i) 



The integrals over q n are conveniently performed successively downwards over Aq n+ i 



q n +i - Qn at fixed q n+1 . The weights \Jg(q n ) = \Jg(q n +i ~ Ag n+i ) require a postpoint ex- 
pansion leading to the naive Jacobian Jq of (\1 14|) and the Jacobian action A e Jo of Eq. (|121|) . 



It goes without saying that the path integral ( |141| ) also has a phase space version. It is 
obtained by omitting all (M/2e)(Ag n ) 2 terms in the short-time actions A e and extending 
the multiple integral by the product of momentum integrals 

dp n 

2nhyJ 'g(q n ) 

When using this expression, all problems which were encountered in the literature with 
canonical transformations of path integrals disappear. 



N+l 

n 

n=l 



( < / R )E„=i 1 k A ?"- e sff'"'(9»W! , «'] 



(142) 
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V. THE PET MODEL IN ONE TIME DIMENSION 



Equipped with thegeneral theory of path integrals in curved spaces we are ready to 
attack the bosonization problem. To become familiar with the subject, consider first a most 
elementary fermion theory described by a Hamiltonian operator 

H= £ -{a'af (143) 

where a\a denote creation and annihilation operators of a fermion at a point. To see the 
difference with respect to boson operators, we shall discuss both options at the same time. 

A. Hilbert Space and Generating Functional 

The states are 

| n ) = _L(at)"|0), n = 0,l,..., (144) 
yn\ 

with energies 

E n = E -n\ (145) 

In the boson case, the quantum number n can run from to infinity, in the fermion case it 
may take only the values and 1, i.e., the energies are 

E = for |0) 

E X = E - for |l) = a t |0). (146) 

The generating functional of all correlation functions of the system is defined by 

Z[r)*, 77] = Tr je-^-* a )f exp i dt{rfa + tfrf) J , (147) 

where T is the time ordering operator and r)(t),r)*(t) are external sources, which are anticom- 
muting Grassmann variables for fermions. The n-point correlation functions are obtained 
from the nth functional derivatives of Z[r)*,rj\. Z[r)*,r)]. 
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The classical Lagrangian of the system is 



L(t) = a*(t)id t a(t) - -\a*{t)a{t)Y 



;i48) 



and the path integral representation for the generating functional ( |147| ) takes the form 



Z[r]*,rj\— / Va*Vaexp 



t a 



dt (L + r]*a + a*r]) 



(149) 



where T is the time ordering operator. For the sake of generality, we first consider a finite 
time interval (t a ,t b ) which will eventually be extended to the entire time axis. The fields 
a*(t),a(t) satisfy periodic of antiperiodic boundary conditions in the bosonic or fermionic 



case: 



a(t b ) = ±a(t a ) 



a*{t b ) = ±a*{t a ) 



(150) 



As long as t b — t a is finite, the generating functional at zero currents Tj(t),rj*(t) is known: 



Z = Z[0,0) = Y,e 



-i(t b -t a )E n 



(151) 



where the summation index runs from n = to infinity for bosons and from to 1 for 
fermions, in accordance with the spectra ( |144j ) and ( |146| ). The expression ( |151| ) is the real- 
time version of the partition function of the system corresponding to an imaginary inverse 
temperature j3 = i(t b — t a ). This follows directly from the spectra (|145|) or ( |146| ), and can 
easily be calculated via path integrals following standard methods (for instance those in 
Chapter 2 in Ref. @). 



B. Collective Quantum Field 



We now introduce a collective quantum field into the path integral via the Hubbard- 
Stratonovich transformation formula [E31 



cxp 



-i f tb dt-[a ] a{t)} 2 
h a 2 



Vp{t) exp \ i I dt 



^-p 2 (t)-p(t)a^a(t) 



(152) 
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which amounts to multiplying (|149|) by the trivial unit factor 

J Vp{t) exp rft - = 1 

with s(t) = so) ait), and integrating out the p-field. Note that because of ( |15U| ), the com- 
posite filed a*(t)a(t), and thus also the field pit) satisfy periodic boundary conditions on the 
interval (t a ,tb). Thus it has the Fourier decomposition 

p{t)=p + p'{t); with p'{t) = Yl (p m e iuJmt + c.c), (153) 

m=±l,±2,... 

with the frequencies u m = 2n/(tb — t a ). The zero-frequency component po is the temporal 
average p = dtp(t)/(t b - t a ); the field p'(t) has a zero average. 

In terms of the Fourier components, the measure path integration for p(t) is 
^ f dp , r dRe p m dim p m 

V P~ I — 11 / , —f . > ( 154 ) 

2ne/iAt m=±i,±2,... J^e/iAt Jire/iAt 



where At = (tj, — t a ). The resulting generating functional Z may be written as 

Z[r)*,r)] = J Va*VaVp 
x exp < i J dt 



a*(t)id t a(t) - p(t)a*(t)a(t) + + r)*(t)a(t) + a*(t)rj(t) 

2e 



, (155) 



where the path integral T>p may be performed by integrating over all Fourier components 
in the standard way. 

Classically, the collective field is proportional to the particle density. Indeed, by extrem- 
izing the action in ( |155| ) we find the relation 

p(t) = ea\t)a(t). (156) 

Integrating out the a*, a fields in ( |155| ) gives 

Z[r)*,r)] = J Vpexp !^A[p) - jf*" dtdt'r)*(t)G p (t,t')r)(t')} (157) 



with the collective field action 

A[p] = ti log Det(G p /i) + f b dt^Q = ±^Tr log^G; 1 ) + f* dt^-, (158) 

Jt a Is H Jt a 2e 
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where G p denotes the Green function of the fundamental particles in an external potential 
p(t), satisfying the differential equation 

[id t -p(t)] G p (t,t') = i5(t-t'). (159) 

This equation may be solved by introducing an auxiliary field 

tplt) = f dt'p'it') + const. (160) 

Jt a 

Inserting the Fourier decomposition ( |153|) we may take 

¥>(*) = E (^e—* + c.c.)= E ^(M^-cc.), (161) 

m=±l,±2,... m=±l,±2,... iU! rn, 

which is a periodic function with a vanishing average. Then we write Eq. (|162|) as 

[idt -Pa- <p(t)\ G p {t, t') = iS{t - t'). (162) 

This is solved by 

G p (t,t') = e-^e^G P0 (t,t'), (163) 

with G Po being the Green function of the fundamental field a(t) for a constant field p(t) = po, 
satisfying the equation 

[id t G po (t,t')-p ]=iS(t-t'), (164) 

and describes the propagation of the fields al o (t) , a po (t) with a Lagrangian 

L Po(t) = al Q it)id t a po {t) - p aJ o (t)a po (t). (165) 

This is the Lagrangian of a harmonic oscillator of frequency u = po- The Green function 
satisfies periodic or antiperiodic boundary conditions in the time interval (t a ,tb) for bosons 
or fermions, respectively. 

For an infinite time interval, the solution of (|164|) is very simple: 

G po (t,t') = e- ip0 ^Q(t - t') (166) 
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for both bosons and fermions. 

For a finite interval, the right-hand side must be made periodic or antiperiodic by adding 
the repetitions, and we find: 

oo 

Gp c) (M') = E {±l) n e~^ t - t '-^ t ^Q{t -t'- (t b - t a )n). (167) 

n=— oo 

The explicit evaluation of the sum on the right-hand side may be restricted to the basic 
interval 

t-t'e [0,t b -t a ), (168) 
where the sum yields in the periodic case 



p -i P0 {t-t') 

G P0 (t,t')= £ e ■ 



,-i P o[t-t'-(t b -t a )n] ' 



2sm[p (t b -t a )/2] 
In the antiperiodic case, we find 



1 _ p—ipa(h—ta) 
n=-oo ^ ° 

e -ipo[t-f-(t b -t a )/2] 

i - . ; —^t, t-t'e [0,t b -t a ). (169) 



-ipo[t-t'-(t b -t a )n] (_\n 



p ipo(t-t') 

G P0 (t,t')= £ e 

n=— oo 
p-i/Oo[t-f-(ti-t )/2] 



t-t' e[o,t b -t a ). (170) 



2cos[p (t fe -t a )/2] 
to be extended outside the interval t G [0,4 — i a ) by antiperiodicity. 

In the original operator language of Eqs. ( |14^) and ( |147| ), the Green function G Po (t,t') 
is equal to the average operator expectation 

G Po (t,t') = (a(t)a\t')) P0 = 1 . — ;/ . (171) 

Tr < e~ ip ° a Tct C*b— *a) k 

For an oscillator state |n), we find an individual quantum mechanical expectation 

n G P0 (t,t') = e -^-*')(n|f (a Po (t)flt o (0) |n> 

= (n + l)6(t-0±ne(* / -*). (172) 

For fermions, only n = and n = 1 contribute. The expectation (171) is obtained by 



averaging these expressions with a pseudo-Boltzmann weig ht factor e -*/w»(*6-**). The result 
coincides, of course, with (|169|) and ( |170| ). 
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The collective field action ( 158|) contains the Tr log of the inverse Green function G p (t, t'). 
To evaluate this, we calculate its functional derivative: 

5 



5p(t) 



±zTr log^G; 1 ) =TGJt,t%, =t+e , 



;i73) 



where the t' —> t limit is specified in such a way that the field p(t) couples to the expectation 



(a\t)a(t)) po = ±(f (a(t)tf(t'))) po \ t , =t+e = ±G p (t,t% l=t+e . 



(174) 



This specification assumes that the terms J^(—dtp(t)a*(t)a(t) + p 2 (t)/2e) in the time-sliced 
version of the path integral (|155|) have the form eJ2 n =i[~ Pn{t)a*(t n )a(t n -i) + p 2 j2e\, with 
the time of a* coming after the time of a (e is the thickness of the time slices). 



For an infinite time interval, the right-hand side of ( 173 ) vanishes trivially due to the 
G-function in ( |166| ). For finite tf, — t a , the right-hand side is nonzero. Inserting the solution 
( |163|) , we see that the (^(independence cancels due to the equality of the time arguments 
and we can replace ( |173| ) by 

6 



6p(t) 



±iTr log^Gr 1 ) = T G P0 (t,t%, =t+e . 



175) 



Due to the constancy of p , the right-hand side is constant. It is equal to the negative 
average particle number n of a harmonic oscillator of frequency po : 

1 



± G po (t,t')\t> =t+e = n= (a f a) P0 = 



Tr [e-^^-^a^a) 



Tr ^e~ ipo ^ ^ tb ~ ta ^ 
Integrating the functional differential equation 

5 



-iP0(tb-ta) 



176) 



5p(t) 



±iTr logger 1 ) 



— n 



(177) 



we find 



± iTr log(iGZ l ) = ±zTr log^'G" 1 ) - n / dtp'{t) 



to 



178) 



The p'(t)-term, however, vanishes due to the periodicity of p'{t), so that ±zTr logliG^ 1 ) 
coincides with ±zTr log(iG~ o 1 ). The associated functional determinant is equal to a real- 
time version of the partition function of a harmonic oscillator of frequency u = po : 
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[Det (G^/i)^ 1 = [Det {G p Ji)\ ±x = Z. 



Pa 




(179) 



This can be written SIS db spectral sum 



-iAipol-l 



y Po 



E- 



(180) 



X _|_ g-iAtpo 

where the summation index n has the same ranges for bosons and fermions as in Eqs. ( |151| ). 
With these results, the generating functional (|157 ) takes the final form 



Z[r)*,r)] 



dpo 



JAtp 2 /2e Z 



2-Ke/iAt 
x / Vip(t) exp 



- / dtip 2 (t) 
Is Jt a 



(181) 



We have changed the integration variables from p'(t) to y?(t). From the measure of p- 
integration (|154|) we see that 

dRe ip m dim <p m 



m=±l,±2,. 



(182) 



since the Fourier components of p'(t) in the integration measure of ( |155| ) and those of ip(t) 
in 



18 1|) are related by p m = iu m (p m . The factors u m are necessary to define the correct 



path integral of a field with a kinetic term <p 2 (t) (see the measure discussion in Ref. 0, 
Section 2.13). Since <p(t) is a massless field, the product of integrals does not include the 
zero-frequency mode of ip(t) — otherwise the partition function would not exist. 
The factors u m are in accordance with the formal functional Jacobian: 



T>p = T>ip det 5(t — t') = const ■ Dip 



183) 



where the constant is the product of all frequency eigenvalues. 

Observe that it is ip(t) which becomes a convenient dynamical plasmon variable, not p(t) 
itself. The original theory has been transformed to a new one involving bosons of zero mass. 
In realistic electron gases they describe plasma excitations JTSI. For this reason, we refer to 
the field <p in the exponent of (|204|) as the plasmon field |15 . 
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VI. COMPARISON BETWEEN ORIGINAL AND BOSONIZED FORMULATIONS 



To see how the bosonization works in detail, let us calculate several properties of the 
model in the two equivalent formulations. 

A. Partition Function 

We begin with the generating functional at zero external currents, the real-time version 
of the quantum partition function. Using the Hamilton operator ( [L43| ) , we have 



Z = Z[0,0] = Tre 



-iAte(ata) 2 /2 



E 



-iAtgn 2 /2 



(184) 



where the summation index has the same ranges for bosons and fermions as in Eqs. (|151|) 
and P0|) . 

The same result is, of course, obtained from the path integral representation ( |149| ) : 



Z = Z[0,0]= / Va*Vaexp 



i / dtL 

. Jt a 



(185) 



if time slicing and measure of integration are defined appropriately 0. 

Consider now the bosonized path integral representation ( 181] ) without external sources, 



Z = Z[0,0] = I d P» e^^ Z f p^exp [f dt*?(t) 
J J2ire/iAt J Vie h a 



(186) 



for bosons and fermions, respectively. Inserting the Fourier representation (|161|) and using 
(|182[) , we see that the path integral over is equal to unity: 



the measure 



T>(f(t) exp 



2e Jt a 



dt(p 2 (t) 



1. 



(187) 



To perform the integral over p , we insert for Z Po the spectral decomposition ( [L80| ), and 



86T) becomes 



Z = Z[0,0] = J 



d P0 p iAtp 2 /2e e ~ iAt P0 t 

2ire/iAt n 



(188) 



After a quadratic completion, the integral over p can be done and yields precisely the 
expression (|184|) . 
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B. Correlation Functions 



For a calculation of the correlation functions of the original fields a*(t) and a(t), we must 
form the functional derivatives of ( |181| ) with respect to the sources rj*(t),r](t), divide the 
result by Z[0, 0], and set the sources equal to zero. Each pair of differentiations 5/5rj*(t) and 
5/5r](t') produces a factor e'^^e^'^Gp^it, t') in the integrand. The path integral over cp- 
fields amounts to calculating the Gaussian averages of these exponentials. For an arbitrary 
functional of <p, these are defined by 



(F[<p]) v = J V V (t)F[ V ]exp 
By Wick's rule, we know that 



i r*b 



2s 



dt(p 2 (t) 



V<p(t) exp 



i r f b 



2e Jt a 



dt(p 2 (t) 



189) 



/ e -M*) e ^(*') 



)<p = ( e 



-i[<p(t)-<p(t')] 



)v = e 2 



\{[v(t)-v(t>)]*) v 



e - \ & 2 (t)> v e - \ iv 2 (t)) Ve (<p(tMt' )) v 



(190) 



where {<p(t)(p(t')) v is the correlation function 



2e 



_ " r -iU!m(t-f) _ b _ 



At 



i\t-t 



'12 



m=l 



- -\t - t'\ + -At. 

2 At 2 1 1 8 



;i9i) 



Hence 



-i<p(i) i<p(t'} 



% (t - t') 2 % 



t-t'\ 



(192) 



2 At 2' 

Note that the t, t'-independent last term in (|191|) has dropped out, so that the correlation 
function of exponentials (e'^^e 1 ^^)^ has a finite limit for At — > oo, in contrast to the 
correlation function of the field ip(t) itself. 

With the result (|192|) it is easy to calculate the correlation function of a boson or a 
fermion field. From (|147|) , its operator expression is given by 



G(t,t') = (fa(t)a\t')) = Z^Tr 



3 -iH(t b -t a ) 



Ta(t)a\t" 



(193) 



Inserting a sum over all intermediate states Sn=o \ n )( n \ = 1; we fi 11 ^ 

oo 

G{t, t') = Z- 1 e- tAtn2/2 e l{t - t,)£n (n + 1), t — t' G [0, t b - t a ). (194) 



n=0 
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The same result is obtained from the bosonic generating functional ( 181|) . For the nor- 
malization factor Z in ( |193| ), this has just been shown. Let us calculate the numerator, 
denoting it by G^{t,t'). Applying to (|181| ) the differentiations 5 2 /5r]*(t)5ri(t'), we obtain its 
path integral 



Gn(M') = / d *> e iAt p ° /2e Z P0 G po (t, t') 



x / V<p(t)e 



-itp{t) i<p{t') 



cxp 



Tel. ^ 



(195) 



The second factor is equal to the correlation function (192). To evaluate the integral over 
p , we write Z po G po (t,t') as a spectral sum 



G P0 ,n(M') = E e 



-iAtp n -ip (t-t') 



(n + 1), t-t' e [0,t b -t a ). 



(196) 



n=0 



After a quadratic completion, the integral over p can be performed and we obtain precisely 
the numerator of ( |194|) of the correlation function. 

For more than one pair of exponential fields e~ lip ^ e t(p ^'\ we have to calculate the ex- 
pectation of functionals of the form exp [i J2% Qi0{ti)] where the numbers qi have the values 
+ 1 for an incoming boson or fermion, and —1 for an outgoing one. The numbers qi may be 
interpreted as the charges of the fundamental fields. After rewriting 



exp 



exp 



dt(p(t)qiS(t - U 



(197) 



we can again apply Wick's rule ( |190| ) and find 

dt<p(t)qi5(t — t, 
1 f°° 

Z J—OD 



(^exp 



(198) 



cxp 



exp 



1 



(199) 



Inserting the correlation function ( |191| ), the right-hand side becomes 

exp -i Qi^j At/16 exp j ~ ^ qiqj[(U - tj) 2 /At - \U - tj 



(200) 
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Since the external sources r)(t),r)*(t) are differentiated pairwise, the total charge q = J2iQi 
vanishes (charge neutrality), so that the first exponential is equal to unity, thus ensuring 
that the expectation has a finite limit for At — » oo: 



cxp 



& 



s i9i ,o exp 



(201) 



It is useful to study the bosonized form of the theory in the operator language to under- 
stand the structure of the Hilbert space. For this it is useful to consider the simpler situation 
of an infinite time interval (corresponding to a zero-temperature equilibrium calculation). 
Then the integral over p i n (|181|) can be done trivially yielding unity and forcing p to be 
zero. The Green function coincides with the vacuum expectation value of the time-ordered 
product 



Go(M') = (0|f (oo(t)oS(O) |o) = e(*-0, 



(202) 



and ( |163| ) yields 



G p (t, t') = e-^V^ett -t'), t> t'. 



(203) 



The generating functional is simply 



Z[r]*,r}] = / Vtp(t)exp 



- / dt^ 2 {t) 



dtdt'tfitfrjit^e-^e^eit - t f ) 



(204) 



To study this theory in the operator language, we take the free plasmon action 

1 fa 



go over to the canonical form 



A 



A=— dtip 2 (t), 
2e ha 



tb dt\p(t)<p(t) - S -p{tf 

ta 2 



(205) 



and identify the Hamiltonian as H = ep 2 /2. After replacing p 
the canonical equal-time commutation rule 

[p(t),m] = -h 
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(206) 

p, ip —> ip, which satisfy 
(207) 



we obtain the Hamilton operator H = ep 2 /2 of the bosonized model. In the Schrodinger 
representation, the operators ip are diagonalized on states \(p) and the functional momentum 
operator p is represented by the differential operator —id /dip. The eigenstates of the Hamil- 
ton operator H consist initially of plane waves which are eigenstates of p with arbitrary real 
eigenvalues p: 

{p\p} = e *<*. (208) 

We are using curly brackets to distinguish the Hilbert space of the <£>-field from that of the 
original a^,a fields. The eigenstates ( |208| ) have the normalization: 

dp{p\p}{p\p'} = 2<K5{p-p'). (209) 



In the operator version, the generating functional (204) reads 



Z ^^ = {bjo} {0|Texp 



dtdt'rf{t)r}{t')e- m) e m,) Q{t - t') 



|0} (210) 



where ip(t) are free field operators. The time-ordered operator on the right-hand side is 
taken between the states of zero-functional momentum. 

We can now generate all Green functions of fundamental particles by forming functional 
derivatives with respect to r]*, 77. First 

5&Z 



(0\Ta(t)tf(t')\0) 



{0|0} 

Inserting the time evolution operator 



5r}*(t)5r)(t') W ^ 

! {0|e-^e^|0}6(t-f). (211) 



e ~iHt = e -iepH/2 ( 212 ) 

the matrix element ( |211|) becomes 

' {o| e ^ ep2/2 e"^ (0) e~ iep2( *"*' )/2 e^ (0) e" i£p2 '' /2 |0} 



{0|0} 



{0 1 e~ im e -^ 2 (*-*')/2 e ^(o) J 0} . (213) 
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But the state e^°)|0} is an eigenstate of p with momentum p — 1, so that (|213|) yields 

— {l|l}e- i£( *-*' )/2 = e-^- 4 ')/ 2 , (214) 
{0|0} 111 ' V ; 



and the Green function (|211 ) becomes 



(0\f a(t)a\t')\0) = e- l£{t - t,)/2 Q{t - t'). (215) 

The same result would, of course, have been obtained for the original fundamental fields 
a^(t),a(t) using the Hamilton operator ( p.43| ): 

(0|f a(t)a + (O|0) = 6(t - t / )(0|e i£{ata)2 * /2 a(0)e- t£(ata)2/2{t -*' ) a t (0)e- i£{ata)2t ' /2 |0) 

= 6(t - t')e- l£{t - t,)/2 . (216) 

Observe that nowhere in the calculation has the Fermi or Bose statistics of the operators 
a(t) and a)(t') been used. This becomes relevant only for higher Green functions. Expanding 
the exponential in (|210|) to the nth order gives 



1 (— V 

Zln] [V*>V] = 77^-+- I dt ^i ■ ■■dt n dt' n r,*(txWt) ■ ■■V*(tnW n 



oo 



{0|0} n l 

x {olTe-^W^ ■ • •e-^ (t " ) e^ (i ™ ) |O}0(t 1 - t[) ■ ■■Q(t n - t' n ). (217) 
The Green function 

(0\fa(h) • ■■a{t n )a\t' n ) ■ ■■a\t' 1 )\0) (218) 
is obtained by forming the derivative 



-i 



,2n 



5 {2n) Z[r]*ri} 



5 V *(t 1 )---S V *(t n )Sr ] (t' n )---5 V (ty 

There are (n!) 2 contributions due to the product rule of differentiation, n\ of them being 
identical thereby canceling the factor 1/n! in fl217|) . The other correspond, from the point of 
view of combinatorics, to all Wick contractions in fl217Q , each contraction being associated 



with a factor (0|e l ^e l ^ t )|0). In addition, the Grassmann nature of source fields r}{t),rf(t) 
causes a minus sign to appear if the contractions deviating by an odd permutation from the 
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natural order 11', 22', 33', .... Denoting a Wick contraction by a common number on top 
of a field operator, we obtain for example 

(OlTa^a^a^a^lO) 

= (0|f a(ti) a(h) l\t' 2 ) a \t[) |0> ± (0|f a(t x ) a(t 2 ) a t(^) a t(^)|o) 



' ■ {0 1 Te - ^* 1 ) e - ^*-* 2 ^ e^' 2 ^ e^ 1 ' 1 ^ 1 } 



{0|0} 

= [9(ta - ti)9(t 2 - t' 2 ) ± e(t a - t 2 )9(t 2 - ti)] (219) 

where the upper sign holds for bosons, the lower for fermions. The lower sign enforces the 
Pauli exclusion principle: If t\ > £ 2 > t' 2 > t[ the two contributions cancel, reflecting the 
fact that no two fermions a? (tt^al (t[) can be created successively on the particle vacuum. 
For bosons one may insert again the time translation operator ( |212|) and complete sets of 
states / <ip|p}{p| = 1 with the result: 

1 



{0|0} 



x {p'\e i ^ i0) e~ i£p " 2/2{t ' 2 ' t ' l) \p''}{p''\e im \0} = e _i£(tl ~* 2)/2 e _i£2( * 2_ *2) e - i£ (*'2-*'i)/ 2 . (220) 



where {0|e _i ^^|p} = 8(1 — p), {p\e~ l ^°^\p'} = S(p + 1 — p') has been used. This again 
agrees with an operator calculation like (|216|) . 

We now understand how the collective quantum field theory works in this model. Its 
Hilbert space consists of states of any functional momenta \p) with p=real. When it comes 
to calculating the Green functions of the fundamental fields of the original theory, however, 
only a small portion of this Hilbert space is used. A fermion can make plasmon transitions 
back and forth between the ground state |0} and the momentum one state |1}, due to 
the anticommutativity of the fermion source fields rj(t),rj*(t). Bosons, on the other hand, 
can connect all states of integer momentum \n}. In either case, the collective-field basis is 
overcomplete as far as the description of the underlying system is concerned. The source 
statistics selects only a small subspace for the dynamics of the fundamental system. 

Note that such a projection is compatible with unitarity. This is guaranteed by the 
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conservation law a^a = const. In higher dimensions, there have to be infinitely many con- 
servation laws (one for every space point) to achieve unitarity. 

VII. NONABELIAN PET MODEL 

We now generalize the above discussion to the nonabelian case and consider a model 
with a classical Lagrangian [compare (|148|) ]. 

1 2 



L{t) = a*(t)id t a(t) - - 



(221) 



and a Hamilton operator 



H 



£ ( ^cr 
a' —a 



(222) 



where a£,, a a with a = 1, 2 denote creation and annihilation operators of a fermion with spin 
up or spin down at a point. 

The generating functional of all correlation functions is 



Z[t)*,t)] = Tr h- iA{tb - ta) f exp 



i / dt(rj*a + a}rf) 



Va*Va exp 



il dt(L + rj*a + a*7]) 

Jt a 



(223) 



in the operator and the path integral formulation, respectively. 

A. The Original Hilbert Space 

To see the difference between fermion and boson systems, we proceed as in the abelian 
case and discuss both options at the same time. The Hamilton operator may be written as 



H 



(224) 



where 



J = a' —a 
2 



(225) 
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is the operator generating spin rotations. These satisfy the commutation rules 



(226) 



The states 



i i\ 

12 5 2/ 



410), = 4|o> 



(227) 



are the basis of a fundamental spin- 1/2 representation of the rotation group. To see the 
transformation properties under finite rotations, we use the fact that every rotation can be 
done with the help of the unitary operator 



U(cp) = e'** 3 . 



(228) 



The right-hand side can be decomposed as follows: 



g-i^ J _ e -iaJ3 e -i/3J2 e -i-yJ3 



(229) 



where ct,/3, 7 are Euler angles. Under a finite rotation, the spin-1/2 operators transform, 
for example, like 



-0j ad t j a= a t os ^ + a t sin ^ 

1 1 2 2 2 

e-^ 2 4e l/3i2 = -4sin^ + 4cos^. 

2 1 2 2 



(230) 



The states have the transformation behavior: 



[/(a, /3, 7) I Ms) =e- W3 e-^e-^|i, S3 ) = |§,4> (e^/V^e"^ 2 ) 



S3S3 



(231) 



where 



d ',.s(^) 



^ cos I — sin I ^ 
V sin f cos f / 



We now form multi-fermion or -boson states 



(232) 
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2s 



n«)|0) (233) 

i=l 

which transform according to higher-spin representations associated with the completely an- 
tisymmetric or symmetric Kronecker products of the fundamental representation (associated 
with all single column- or row-like Young tableaux). A system with two spin | particles has 
spin for fermions and spin one for bosons. Three-particle states vanish for fermions and 
have spin 3/2 for bosons. In the bosonic case, 2s spin 1/2 particles couple to spin s. 

Explicitly, the properly normalized states of total spin s and magnetic quantum number 
m are given by 

\s,m) = = (a\r m (atr m \0)- (234) 

J (s — m)\(s + my. 

Under finite rotations e~ tlf ' 3 , they transform like 

e-^\jm)= £ |jm')i^, m (a,/3,7)= £ \jm')e^™'^(jm'\e^\jm), (235) 

m'=—j m'=—j 

where 

dL'JP) = (jm'\e- iJ ^\jm) (236) 

is given by 



4n,M = 



(j + m')\(j - m')\ ^ 
\ (j + m)\(j-m)\ ^ 



j + m 



j — m 
k , 



yj -rri -k J 

2k+m'+m / 2j—2k—m'—m 



(n\ 4K-rin Tin / n 
cos|J ^sin|j (237) 

From the above analysis it is obvious that the real-time partition function of the model 
has the spectral sum 

Z = Z[0, 0] = E( 2 J + l)e~ £jU+1)/2 . (238) 
j 

In the bosonic case, each spin j = 0, ±1/2, ±1, . . . occurs precisely once with (2j + 1) 
orientations m = —j, . . . , j. In the fermionic case, only the spins j = 0, ±1/2 occur. 
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B. Collective Quantum Field 



Let us now bosonize the theory ( |223| ). A collective vector quantum field p is introduced 
into the path integral representation ( [223| ) via a Hubbard-Stratonovich formula analogous 
to (PI: 

1 



r°° , e 
exp ^ —i I dt- 

-oo 2 



a*-a(t) 



y vp(t) exp |z y"°° 



2s p 2 (t) - p(t)a* -a(t) 



(239) 



which amounts to multiplying 



by the trivial unit factor 



Vp(t)exv\— dt[p{t)-v{t)f 



2e Jt a 

with v(t) = eaf cra(t) /2, and integrating out the p-field. For an infinite time interval At, 
the integral over the temporal average p = J^ b p(t)/(t b — t a ) of the collective field is forced 
to be zero as in the abelian path integral (|181|) . Then the generating functional is simply 



Z[i]\i]\ = J Dp' exp ju[p'] - J°° dtdt'r)*{t)G p ,{t,t')r){t' 



(240) 



where p'{t) has no temporal average and the Green function G p >(t, t') satisfies the differential 
equation 



[id t - p'(t) ■ a/2] G p/ (t,t') = %5{t - t'). 



(241) 



This equation may be solved by introducing an auxiliary 2x2 hermitian matrix field $(t) 
(p ■ cr/2 via the following identity 



e -i*(t) = f e -ift eo & Si*)-'/* 



(242) 



in terms of which 



G p ,{t,t') = e-^Goit - t'y*^ = e -^*W e »(*')e( f _ if) 



(243) 



thus generalizing ( p.63| ) and ( |203| ). 

We now calculate the Tr log term in ( |240| ) . From ( |241| ) we see that 
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5 1 

r ±zTr \og(iG p ^)] = T -tr [<TiG„(t, t')]\ v=t+e = 



5p'(t) 



(244) 



where the t' —> t limit is specified as in the abelian case [see ( |174|) 1. Inserting the solution 



( j243|) , we find we see that the B-function in (|203|) makes the functional derivative vanish 
and the Tr log becomes an irrelevant constant. 

Note that for a finite time interval (t a ,t(,), the functional properties of abelian and non- 
abelian models are quite different from each other. Then ( |244j ) becomes 



8 



5p{t) 



±*Tr log^ 1 )] = ^tr[^-WG Po (M')]|^+e 



0. 



(245) 



Due to the presence of the er-matrix, the Euler angles do not disappear from the right-hand 



side, in contrast to ( p. 771 ) P7 



Returning to the case of an infinite time interval (t a ,tb), the generating functional is 

\ 1 2 



Z\n .//] • / Dp(t)r W {^-J^dttT 



£ e -i*(t) \ e »(t) 

dt ) 



d^*(t)r/(t , )e" 4 * W e i * ( '' ) e(t - t') 



(246) 



At this place, we observe another important difference with respect to the abelian case. 
There, the kinetic term in the exponent could simply be rewritten as <f 2 {t). Here, this is no 
longer possible. The kinetic term contains interactions between the three field components. 
In order to exhibit these in a familiar form, we express e 1 ®^ in terms of Euler angles. This 
defines the 2x2 unitary matrix 



-i#(t) _ 



u{a(t),p(t)n(t)). 



(247) 



The kinetic term in the action ( |246| ) can then be rewritten as 



tr 



dt 



e -i«(t) e »(t) 



tviuu-^aitimnmf 



Inserting for U(a(t), /3(t), 7(t)) the explicit Euler angle form as in (|231| ), 



(248) 



U(a(t),P(t),7(t)) = e -^3/2 e - i/3ff2/2 e- i7CT3/2 , 



(249) 
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we find that the three components of p(t) coincide with the components of the angular 
velocities of a spinning top whose orientation is described by the Euler angles a, f3, 7: 

Pi(t) = ujiit) = —f3 sin 7 + a sin (3 cos 7, 
p 2 {t) = oo 2 (t) = $ COS7 + asm/3 sinj, 

p 3 (t) =u 3 (t) = a cos (3 + j. (250) 

The generating functional can therefore be rewritten in terms of Euler angles as follows: 

Z[ri*,r]}= I VaV cos f3V-fFexp (i- [°° dt u 2 (t) 
J I 2e J —00 

dtdt'r]*(t)U(a{t), Pit), 7 (t))C/t( a (t'), P(t'), ^{t')) V {t')Q{t - t')\ . (251) 



Here F is a functional Jacobian arising when changing the integration variables p(t) to the 
invariant measure in the space of Euler angles at(t), /3(t), "f(t). 

C. Measure of Integration in Bosonized Theory 

At this point, the new results on variable changes in path integrals in Ref. come into 
play. These variable changes are governed by the quantum equivalence principle. Let us first 
introduce a trivial change of integration variables from p(t) to variables 

t 

1 



00 



Q(t) = I dt'p(t'). (252) 

We can then rewrite / T>p(t) as 

' VQ(t). (253) 



In Eq. ( |250| ) we have seen that pi(t) coincide with the components Ui(t) of the angular 
velocity. These are linear combinations of the Lagrangian velocities q^(t) = (6c(t), $(t), j(t)). 
There exists the following relation between the velocities Q*(t) and q^(t): 

with the matrix 
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1 ■ R ■ n \ 

sin p cos 7 — sin 7 (J 



V(a(*), /?(*), 7(*)) 



Equation 



sin /3 sin 7 cos 7 I • (255) 
cos /3 1 

is a nonholonomic mapping of all paths in the parameter space of Euler 



V 



angles into paths Q(t). The former space has a constant curvature, the latter space has 
no curvature, but a nonzero torsion [14, j|. For a finite time interval (t a ,tb), the mapping 
follows the integral equation (|67|): 



= q"(ta) + / dt'efiqit'))®^'). 



ta 



(256) 



(Qt|Q't' 



iM(AQ) 2 /2e 



(257) 



According to Eq. ( |129|) , the correct path integral in a space with curvature and torsion 
is found as follows: In a flat-space with cartesian coordinates Q, the path integral is known 
to have the time-sliced form: 

1 N r roc 1 N+l 

D U / d D AQ n I]« 
<J2meh/M »=i L 00 J »=i 

where the coordinate differences AQ n = Q n — Q n _i appear in the exponent and in the 
time-sliced measure. This measure corresponds directly to the naive time-sliced version of 
the measure ( |253| ) in the present model. 

N+l 



f vQ(t) -> n d ° A Q 

J n=2 



(258) 



The coordinate differences AQ l n are now mapped into a space with curvature and torsion via 
the nonholonomic mapping (|256|) , which is uniquely carried out along the classical short-time 
trajectories. Under this mapping, the short-time actions go over into the actions calculated 
along the classical trajectories, just as postulated in curved spaces by DeWitt |jj (who 
followed in this respect the original observation by Dirac 0, from which Feynman derived 
his path integral representation). As emphasized above, the classical trajectories in the 



presence of torsion are autoparallels, not geodesies [13 



The image of the path measure in g-space is according to (O 
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1 



N 

'2mhe/M n=i 
with an effective potential 



/ d D q n 



9{Qn) 



-D 



2meh/M 



x exp 



„• AT+l 

- £ (A e + eV ea ) 

11 n=l 



(259) 



6M 



(260) 



where the curvature scalar R is defined by the contraction R = g vX R u \ of the Ricci tensor. 

Inserting the Euler angles for q 11 , we may write the measure in the generating functional 
Q as 

VaVcosPV-ye£« dtVefl . (261) 



The action is time-sliced as follows: According to Ref. 0, Section 8.10, one first defines 
a sliced action near the spinning top 



A 



N 



- N+l 

C 71=1 



1 - -tT{U n U-\) 



with 



(262) 



(263) 



The path integral Q251|) without the external currents can then be solved exactly. The action 



fl262Q is not yet the correct one, due to the fact that the differences in (|262|) do not measure 
the sliced geodesic distances. A geodesic correction must be applied which is of fourth order 
in Ag M , as explained in Ref. 0, Section 8.9. 

After this, we calculate (see Ref. 0, Section 8.11) 

Z[0,0] = lim V(2j + if e -i(t b -t a )ejU+i)/2 = L ( 2g4 ) 

t b -t a —>00 Z — ' 

There is no extra term proportional to i? as in DeWitt's path integral for the spinning top. 
It is the quantum partition function of a spinning top in the limit t —>■ oo, where only the 
ground state survives. Note that there is no extra term proportional to R as in DeWitt's 
path integral for the spinning top. 

If we add the external currents, each derivative with respect to r]*(t) or rj(t') produces a 
factor U(a(t), fi(t), j{t)) or [/ _1 (a(t), (3(t), j{t)) in the integrand, respectively. 
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VIII. HILBERT SPACE OF BOSONIZED NONABELIAN MODEL 



In the abelian case, the Green functions of the initial bosons or fermions did not involve 
the full Hilbert space of the bosonized theory. The same thing is true in the nonabelian 
case. The initial particles are represented only by a subset of the wave functions of the 
spinning top. This is seen by calculating the two-point correlation function, obtained from 
the functional derivatives S 2 / 8rf(t)8r](t') of the generating functional Z[q*,r]]. 

In the operator form (|223|) of the generating functional, the two-point correlation function 
is given by the expectation value 

Gw(M') = (0\fa m (t)al f (t')\0), (265) 

for which we easily calculate 

G mm ,{t,t') = 5 mm ,e- iAE ^Q(t - t'), (266) 

where AE is the energy difference between a state carrying one boson or fermion and the 
vacuum state |0): 

AE = 3e/8. (267) 

In the bosonized theory we differentiate ( |251| ) and find 
Gw(M') = / VaVcMpD-yj^w* eX p \± f°° dtu, 2 (t)} [U(t)W(t)] mm ,Q(t - t'), 

J L J — oo J 

(268) 

with U{t) short for U(a(t), f3(t), j(t)). 

As in the abelian case, we evaluate the bosonized expression ( |268|) in the operator lan- 
guage using the Schrodinger representation. Due to the presence of the correction factor 
e ijdtv ctl m measure of the path integral ( |268| ), the Hamilton operator associated with 



the action in ( |268| ) is proportional to the Laplace-Beltrami operator 



A = -Ld^g^, (269) 
V 9 
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where g^ u is the inverse of the metric defined by the kinetic term in the classical La- 
grangian having the form 



In our model 



2e 



1 1 cos/3^ 



X 



(270) 



/ 



1 

cos/3 1 

The Hamilton operator contains no extra term proportional to the curvature scalar, and 
coincides with the one arising from quantizing the generators of the rotation group in the 
classical expression 



H = -J 2 



(271) 



leading to the well-known operator 



H 



dp + cot 



Pdp+ (l + cot 2 /?)<9 7 



2 _ 2cos^ 



sin 2 p a sin 2 /3 



(272) 



This was shown in Ref. 0. 
The eigenfunctions are 



{a(3~{\jmm'} = D 3 mm ,(a, ft, 7), 



(273) 



with the energies 



Ejmm' = 77 j (j + 1) 



(274) 



In this Schrodinger representation, the correlation function ( |268| ) is given by the expec- 
tation value 



Gw(M') = {0|D^(*)D 1 / 2 *(Om'fc|0}e(t-f) 



(275) 



where we have replaced the matrices U(a(t), f3(t), j(t)) by the spin- 1/2 representation ma- 
trices D^,(at(t), P(t), j(t)) of Eq. (|235|) , and written them short as D^, (t), as we did with 
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U(t). The vacuum state has the Schrodinger wave function {a, (3, j\0} = -D$ (a,/3, 7) 



1 / V8n 2 , and an energy 

£0,0,0 = 0. (276) 
Inserting the time evolution operator, we write 

D(a(t),f3(t),l(t)) = e i6t D(a(0),m,im^ iM (277) 
with H of (|272|) and find a phase 

e -iA£(t-t') ) (278) 

where AE is the energy difference between the boson wave function \ jmm'} = 1 1/2, 1/2, 1/2} 
and the ground state |0} = |0,0,0}. Its value is the same as in the operator calculation 

Then ( f275| ) reduces to the integral 

Grnmi(t, t') = ^2 docd cos /3c?7 
k J 

x {0|a/M^W,7)S^^ (279) 
Using the unitarity property of the rotation functions Z) 1//2 (a, /3, 7) 

DlHia, (3, 7 )^>, /3, 7) = <W, (280) 

we can rewrite this as 

G mm ,(t,t') = 8 mm > J dadcos^d7{0|a^7}{a^7|0}e- iAf;(t -*' ) e(t -t') 

= 6 mm ,e- iAE We(t-t'), (281) 

which is, of course, the same as in (|266|) . 

In this expression we observe a nonabelian version of the projective properties of the 
bosonized theory in the Hilbert space of all rotational wave functions. At the level of spin 
1/2, there are four rotational wave functions D±y 2 ±1 ^ 2 (a, (3, 7). The correlation function 
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( |281|) , however, contains one contracted index which makes the angle 7 disappear. The 
same happens in all higher-point correlation functions. Thus, the correlation functions of 
the bosonized theory make use only of a subspace of the total Hilbert space of the spinning 



top in which the Euler angle 7 is absent. The correlation function ( |281| ) looks as though 

1 /2 

the wave function of a spin-1/2 particle were ^(a,/3, 7) oc ^2 k D k ' ±1 , 2 {a, f3,j). These are 
orthogonal and complete in the scalar product defined by 



27T /*7T /"27T 

Jo 



8tt 2 



' d0id P Sil1 P d 1 D m[mi ( tt > & i) D Z' 2 ra 2 («> & 7) = ^m' x m' 2 S mim2 6 jlj2 — - - . (282) 

This subspace of top wave functions is equivalent to the space of spherical harmonics 



Yi m (j3,a) = \j2l + l)/47rZ)^ l0 (a, (3, 7). Except for the presence of half-integer spins, the 
spectrum corresponds to that of a particle on the surface of a three-dimensional sphere, 
where the energy eigenvalues ej(j + l)/2 appear only (2j + l)-times rather than (2j + l) 2 - 
times in the spinning top. This is the selection mechanism reducing the partition function 
of the spinning top ( |264| ) to the smaller sum ( |238| ) over the initial states. 

If the initial fundamental particles are fermions, the orthogonality relation of the rota- 
tion functions _D 1//2 (a, (3, 7) together with the Grassmann algebra ensure that the bosonized 
theory represents properly the anticommutation rules of the original fermion operators. 

If one wants bosonized particles to cover a Hilbert space that is completely equivalent 
to the spinning top, one must start with twice as many bosons as before. The appropriate 
Lagrangian is then 



L{t) = a*(t)id t a(t) + b*(t)id t b(t) - - 



a*(t)^a(t) + b*(t)^b(t) 



(283) 



and the Hamilton operator 



H 



a' —a 



+ b^b 



(284) 



This can be written as 



2 



(285) 
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where 



j(i)^ a t-a, j( 2 > ^ &t^s 

2 ' 2 



(286) 



are two independent sets of angular momentum operators with the commutation rules 



J 1 J 2 



J 1 J 1 

J l 7 J 



J 2 J 2 



ufe^fe 7 



(287) 



The Hilbert space consists of the states 
\nlnln\n b 2 ) = ' 



n ( \\n%\ri}\ \n\\ 



(4)^(4)^(61)^(^)^10). 



(288) 



H : "'2 : "'l : "'2 : 

If we consider only the states with an equal number of a and b particles, 



(a f a- 6+6)1^) = 0, 



(289) 



the Hilbert space is equivalent to that of the spinning top. To enforce ( |289| ), we have to 
extend the Lagrangian ( |283| ) by a Lagrange multiplier 



X(t)[a*(t)a(t) - b*(t)b(t)]. 



(290) 



It is worth pointing out, that a free-oscillator version of the Lagrangian ( |283| ) with the 
constraint fl290|) . 



L(t) = a*{t)id t a{t) + b*{t)id t b{t) - uo [a*{t)a{t) + 6*(*)6(*)] + \{t)[a*{t)a{t) - b*(t)b(t)}, (291) 

arises from a nonholonomic transformation of the path integral of the hydrogen atom (see 
Chapter 13 in 0). Thus, the path integral of the hydrogen atom could, in principle, also 
be solved by a Duru-Kleinert transformation to that of a spinning top containing an extra 
energy term proportional to a*(t)a(t) + b*(t)b(t). 
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IX. NONABELIAN VERSION OF HUBBARD-STRATONOVICH 
TRANSFORMATION FORMULA 



A crucial role in the bosonization procedure is played by the Hubbard- Stratonovich 



transformation ( [239| ). After replacing p by Q according to ( 252 ) and performing the non- 
holonomic transformation Q254J) to the Euler angles, this can be rewritten as 



exp < — % 



OO £ 

dt- 

-OO 2 



a\t)-a{t) 



T>otD cos (3T>^j e 

1 



i S- oa dtV 'B 



x exp J dt 
Equivalently, there exists the following nonabelian identity: 

' VaV cos/3P 7 e^-°°- dtV * e exp [°° dt [w(t) - v(t)] 2 



, l£ " \t)-u{t)a\t)-a{t) 



(292) 



(293) 



valid for an arbitrary time-dependent vector field v(£). The time slicing of the action has 
to be done as in Eq. ( |262|) with the subsequent geodesic correction explained in Ref. 0, 
Section 8.9. 

For a finite time interval (t a ,tb) these formulas contain, of course, an extra integration 
over the zero mode of the initial collective quantum field p(t), as in ([[]): 



dp 



,iAtp 2 Q /2e 



'2TTE/iAt 

The proof of formula ( |293| ) is quite simple: We take any time-dependent matrix U v (t) 
solving the differential equation 



U v (t)U~\t) 



-v • a. 

2 



and rewrite the exponent in ( |293| ) as 

' d 



U i-[U v (t)U(t)} [UvWit)]- 1 



(294) 



(295) 



Changing variables from the Euler angles of U(t) to those of U v (t)U(t), and using the 
invariance of the integration measure under this group operation, we obtain directly the 
independence of the path integral (|293|) of v(£). The normalization to unit is trivial. 



Generalizations of this formula should be useful in bosonizing other nonabelian theories. 
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X. CONCLUSION 



The bosonization of the simple spin model requires taking proper care of the nontrivial 
Jacobian which arises by the nonholonomic field transformation to the Euler angles. Thus, 
in addition to the solution of the path integral of the hydrogen atom, bosonization is a 
second important example for the power of nonholonomic field transformations in relating 
path integrals of completely different systems to each other. The nontrivial Jacobian arising 
in the transformation process is uniquely derived from the quantum equivalence principle. 
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FIGURE HEADINGS 

Fig. 1: Crystal with dislocation and disclination generated by nonholonomic coordinate 
transformations from an ideal crystal. Geometrically, the former transformation introduces 
torsion, the latter curvature. 

Fig. 2: Images under a holonomic and a nonholonomic mapping of a fundamental path 
variation. In the holonomic case, the paths x(t) and x(t) + Sx(t) in (a) turn into the paths 
q(t) and q(t) + Sq(t) in (b). In the nonholonomic case with S^ x ^ 0, they go over into q(t) 
and q(t) + <5q(t) shown in (c) with a closure failure 6 M at tf> analogous to the Burgers vector 
6 M in a solid with dislocations. 
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